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Abstract 



We consider the quantum inverse scattering method for several mixed integrable mod- 
els based on the rational SU{N) i?-matrix with general toroidal boundary conditions. 
This includes systems whose Hilbert spaces are invariant by the discrete representations 
of the group SU{2) and the non-compact group SU{1, 1) as well as the conjugate repre- 
sentation of the SU{N) symmetry. Introducing certain transformations on the quantum 
spaces we are able to solve generalized impurity problems including those related to 
singular matrices. 
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1 Introduction 

The quantum version of the inverse scattering method is nowadays accepted as the foundation 
of the algebraic theory of integrable models of quantum field theory and classical statistical 
mechanics in (1 + 1) dimensions P3I21- This approach has made possible not only the discovery 
of new important models such as the Heisenberg chain with arbitrary spin [HI E] but also 
precipitated the notion of quantum group symmetry |S1IE]- 

An essential ingredient of the quantum inverse scattering method is frequently denominated 
the Yang-Baxter algebra whose generators can be thought as the elements of square N x N 
matrix T^{X). The symbol A emphasizes such A^- dimensional auxiliary space and A denotes 
the spectral parameter. This associative algebra is then generated by the following quadratic 
relations 

^(A - i2)TAX) ® T^(/x) = T_^{fi) ® X4(A)i?(A - /i), (1) 

where -R(A) is a A^^ x A^^ matrix whose elements are complex numbers. Here we are assuming 
that this i?-matrix is additive with respect the spectral parameters. 

The i?-matrix -R(A) is required to satisfy a sufficient condition, that guarantees the asso- 
ciativity of the algebra (^, known as the Yang-Baxter equation 

i?12(A)/?23(A + /i)i?12(/^) = i?23(/i)^12(A + /i)/223(A), (2) 

where RabW denotes the action of the i?-matrix on the tensor product space Aa (8 Ab- 
It is not difficult to see that the trace of T4(A) over the auxihary space 

T(A) = Tr^ [T^X]] , (3) 

gives origin to a commutative family of operators, i.e. [T(A),T(/i)] = for arbitrary values of 
A and /i. In this sense, T(A) can be regarded as the generating function of quantum integrals 
of motion. 

One possible representation of T^^iX) is provided by numeric N xN matrix Qj^ whose entries 
do not depend on A and that satisfies the relation 

[R{x),gA^GA\ =0, (4) 
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and therefore being direct related to the underlying symmetries of R{X). 

General representations of (Q) for a given i?-matrix depend on the spectral parameter A 
and are often called /^-operators C^j{X). They are viewed as matrices on the auxiliary space 
A whose elements are operators on another space Vj known as quantum space. The simplest 
one occurs when the auxiliary space A and the quantum space Vj are isomorphic, since in this 
case Eq.Q becomes equivalent to Eq.((T} provided we set 

jCAjW = PAjm, (5) 

where P^j is the exchange operator on the space A^Vj. 

An extremely important property of the Yang-Baxter algebra is that the tensor product of 

L 

two representations is still a representation of this algebra. For example, let C^j ^ TT ®^' for 

L distinct /^-operators which satisfy the Yang-Baxter algebra with the same i?-matrix. Then 
the following ordered tensor product of /^-operators, denominated monodromy matrix PI Ej 

TaW = Ga^alWCal-iW ■ --CAiiX), (6) 

is also a representation of the same Yang-Baxter algebra. 

In the context of classical vertex models of statistical mechanics Caj{X) represents the 
possible Boltzmann weights at the local site j of a square lattice of size L, Qa V^^Y the role 
of possible toroidal boundary conditions compatible with integrability [3 [S] and the operator 
T(A) turns out to be the corresponding row-to-row transfer matrix. The purpose of this 
paper is to investigate the eigenvalue problem for the transfer matrix T(A) of such mixed 
vertex models whose underlying i?-matrix is the simplest rational solution of the Yang-Baxter 
equation invariant by the fundamental representation of SU{N), namely PlITU] 

N N 

Ra,{x) = V E ^- ® ^- + ^ E ^S ® 4- (7) 

o=l «,/3=l 

where e)^^ are the N x N Weyl matrices acting on the space Aa and rj is the so-called quasi- 
classical parameter. It turns out that the admissible boundary representations Qa for this 



i?-niatrix are arbitrary N x N matrices that will be represented by 

N 

Ga= Yl 9ape^fJ- (8) 

o,/3=l 

The purpose of this paper is twofold. On the one hand, we extend our recent efforts 
[TT] in solving fundamental SU{N) vertex models with rather general non-diagonal toroidal 
boundary conditions to other interesting group invariance of the quantum space that are not 
isomorphic to the fundamental representation of the SU{N) symmetry. Secondly, we explore 
the consequences of the existence of transformations on the quantum space to clarify possible 
relations between the eigenvectors of T(A) with diagonal and non-diagonal boundary condi- 
tions. This fact makes possible the solution of generalized mixed vertex models that combine 
both boundary and spectral dependent representations at any site of the chain. 

This work has been organized as follows. In section |21 we present the algebraic Bethe ansatz 
solution of SU{2) mixed vertex models when their Hilbert spaces are invariant by higher spin 
representations of SU{2) and the discrete D^{k) representation of the non-compact group 
SU{1, 1). In sectioniniwe discuss similar analysis for the isotropic vertex model that mixes the 
fundamental and the conjugate representations of SU{N). In section |3] we show that suitable 
transformations on the quantum space are essential for the solution of general classes of mixed 
vertex models. This includes those whose boundary representations are singular as well as the 
integrable spin-S* Heisenberg model with the most arbitrary non-diagonal boundary conditions. 
The appendix A is reserved for the study of the completeness of the Hilbert space of some of 
the mixed vertex models mentioned above with L = 2. In appendix B we present the explicit 
expressions of the quantum transformations for all the mixed systems described in this paper. 

2 Mixed SU{2) vertex models 

The purpose of this section is to solve the eigenvalues problem for the transfer matrix (jS)) of 
mixed vertex models whose underlying i?-matrix ((7j) is SU{2) invariant. One natural way of 
producing such mixed vertex model is by choosing /^-operators intertwining between general 



representations of SU{2). The other manner is to look for reahzations of the Yang-Baxter 
algebra (^ in terms of different algebraic structures for the quantum space Vj, even those 
based on non-compact groups. In what follows we shall explore these both possibilities. 

2.1 The higher spin reaUzation 

The approach of using arbitrary representations of SU{2) to built up mixed vertex models goes 
probably back to the work by Andrei and Johannesson ^2] who studied the spin-1/2 Heisen- 
berg model in the presence of an impurity of spin-S". Latter on, de Vega and Woyanorovich JH] 
have used similar idea to construct integrable Heisenberg chains with alternating spins resem- 
bling ferrimagnetic models. However, those works have assumed explicitly periodic boundary 
conditions and therefore not dealt with the most general forms of integrable impurities that 
are going to be considered here. 

The higher spin realization of (P) is built from the local /^-operators |31 E] 

\ 'nsj [>- + i]-r,Si, ) 

which act in the tensor product of the auxiliary space of spin-1/2 and the space Vj carrying 
the spin-S* representation at the j-th site. The construction of the transfer matrix is standard 
as described in the introduction and it is given by 

t[%-'^\\) = Tr^ [gAJO.±''\x) . . . C'jf^\x)] . (10) 

In order to diagonalize this operator by the quantum inverse scattering method it is desirable 
the existence of local vacuum vectors |0) . such that the action of the monodromy matrix Tj^iX) 
on the global state 

L 

io)=n®io).' (11) 

gives as result a triangular matrix. 

Since the boundary matrix ^^ is arbitrary, the tensor product of the standard highest spin- 
Sj states does not work as an appropriate global reference state. Following our previous work 



[TT] , suitable local states can be found by introducing a set of Baxter's gauge transformations 



M, 



on the spin-1/2 auxiliary space such that 



C%'^\\) = M-l,C%'^\\)M,, 



(12) 



where Mj is an invertible 2x2 matrix whose entries are denoted by 



M, 



^j ^j 



Vj ^i 



(13) 



In terms of these new /^-operators transfer matrix ()1()|1 becomes 



rp{Sl,...,SL) I 

^1/2 



(A) = Tr^ M^'G^Ml+iT^X) 



(14) 



where the Xa{X) is the transformed monodromy matrix T^^iX) = CJ^ ^ {X) ■ ■ ■ C^i ^ (A). 

~(- s) 
We now search for gauge transformations Mj in such way that £_P' " (A) is annihilated by 

say its lower left element for general values of the spectral parameter A. This property turns 

out to be equivalent to the following condition for the ratio n,- = — 



>-l 



{pj+1 - Pj) - vS]{Pj+i + Pj) + vSjPj+iPj - vSj \ |o); 



ASj) 



0. 



(15) 



(s,) ■ 



This problem can be solved by expanding the vector |0)^- in terms of simultaneous eigen- 



kets \Sj, rrij) of 5*1 and 5*1 



|0)f)= ^ C7(5„m,)|^,,m,). 



(16) 



mi=—Si 



Substituting this ansatz in Eq.()15p and by taking into account the well known action of 
the angular momenta ladder operators S^ on the basis \Sj,mj) we find (25* + 1) restrictions 
for the coefficients C{Sj,mj) given by 



^ [-(Pj+i + Pj)mjC{Sj, rrij) + PjPj+iC{Sj, rrij - 1) J{Sj + mj){Sj - m^ + 1 



'mi=—Si 



(^ + 2)(Pi+i " Pj)'^iSj^ ^j) - ^'(•^i' ^J + 1) V (•S'i - mj){Sj + nij + 1)] 15^, mj) = 0. (17) 



A necessary condition for the solution of these equations such that |0)^- is independent 
of the spectral parameter occurs when pj+i = pj and therefore the ratio — for j = 1, . . . , L is 
keep fixed. If we substitute this condition in Eq.(|17|) we end up with recurrence relations for 
the coefficients C{Sj, rrij) which can be solved in terms of an arbitrary normalization C{Sj, Sj) 
constant. The final result for the local vacuum |0)^- can then be written as follows 



io>f'-^E^^ „,,4),(^_,„„ ^r^^(^.^)i^.--.) 



and the action of the operator C_^j (A) in this state is given by 



{>^ + 'i-vS,)'-^^^\ot' 



^\>^m'= ' ' ' I'y^^^^'^ ^^ _ _ - ^^^^^^^ ^^^^^^^ 1 (19) 



where the symbol # denotes non-null states. 

The constant ratio pj is then selected out by imposing that the transformed boundary 
M^^Q_aMl+i matrix becomes a diagonal matrix. This restriction lead us to two possible 
values for pj = p^'^\ namely 

P^ ' = ^ , (20) 

and consequently with the help of Eqs. ()lH18p to two choices for the global reference state 

|o>w. 

We now have the basic ingredients to turn to the algebraic Bethe ansatz solution of the 
eigenvalue problem associated to T^,^''"' ^ (A). The first step is to seek for a convenient repre- 
sentation for the monodromy matrix T^{X). Previous experience in dealing with the quantum 
inverse scattering method for two-dimensional auxiliary spaces [H |2] suggest us to adopt the 
form 

ZW = I t^^^ ^J^^ I • (21) 

\C{X) D{X) ) 

The property (jl9|) help us to identify the elements of the transformed monodromy matrix 

that acts as particle creation and annihilation fields over the pseudo- vacuum |0)^ ' . We see 



that -B(A) play the role of creation field and C{X) is an annihilation operator thanks to the 
property C'(A) |0) = 0. On the other hand, the diagonal fields satisfy the following important 
relations 

^(A)|0)(^) = fl(x + UvS,)-^\0)^^\ (22) 



L 



5(A) |0)<- ^ n(^^i-.-.)^;,^^|0)-. (23) 

(±) 



From these properties and the help of Eq.(j2ni) it is not difficult to see that the states |0) 

/ 



are themselves an eigenstate of T^,2''' ^ (A) with eigenvalue Aq^.^''"' ^ (A) given by 



L 



a£-'^^)(A) = g[% n (^ + i + ^^.) + 4^ n (^ + i - ^^.■) ' (24) 

where the phase factors g[2 are the eigenvalues of the matrix ^_4 

(J 

^1/2 - 2 



(±) (fi-ll + 922) ± \/{gil - 922^ + ^912921 .„p.s 

9\/2 = ^ • (25) 



This result strongly suggest us that the arbitrary eigenvectors 10)^ ' oiT^,^''"' ^ (A) can be 
put in the form 

n± 

|0)W=n^(Af)|O)(^). (26) 

The symmetry [-R(A), Mj ® Mj] = implies that the gauge transformed monodromy matrix 
Xa{X) satisfies the Yang-Baxter algebra (0) with the same i?-matrix as the original monodromy 
matrix T^{X). As consequence of that, the matrix elements of T^{X) satisfies the same set of 
commutation rules of the periodic six-vertex model PQ Ej and from now on the main steps in 
the eigenvalue solution of T^,^'"' ^ W become very similar to that of this well known vertex 
system. Considering that these details have appeared in many different places in the literature, 
see for instance ref. P3I2], here we shall present only the final results. By imposing that |0) 
are eigenstates of Ty^''"' ^ (A) and performing the convenient displacement A^- ^ ^ A^- — | 



we find that the corresponding eigenvalues A„j!''' ^ (A) are given by the expression 




j=l i=l \ ^2 

where the rapidities A^- satisfy the following system of transcendental equations 

yi/2 T-T \ - \ +v , . 

-J±)1K(±) .(±) ^- (28) 

fi'i/2 ;=i \ ~ \ ~V 

We would like to close this section with the following comments. First we emphasize that 

L 

the integers n± play the role of particle numbers sectors satisfying the constraint n± <2\ Si. 

i=l 

In appendix A, we illustrate this fact by presenting the details of a study of the completeness of 
the eigenspectrum ()27p28|) for L = 2. Our final observation concerns with physically interesting 

spin chains that commutes with the transfer matrix ()1()|1 . Though the integrability does not 

(- s) 
depend on how we distribute the ^-operators CJj ^ (A), the construction of local conserved 

charges commuting with T^^''"' ^ (A) does. One interesting case is when we have only one 

impurity i2-operator C_^l (A) sitting at the end of the chain of L — 1 spin-1/2 C^^^ (A) 

(- -) 
Boltzmann weights. Because -C A'^ (A) is proportional to the exchange operator P_aj at the 



Tljr"'\^) 



-1/2 



around this point. 



point A = 0, we can produce local charges by expanding In 

The second term in this expansion is the associated Hamiltonian and its general expression 

reads 

L— 2 

\ ^ I r^2'2) /'n\\ I I r^2'^) /'n\\ i ( r^2'^> r r\\\ r>-l I r^2'2> /n\ \ r> r^2'^) 



^ = -E (4^(0))" +(^K!l(0))' +- {cti]My Ql-i (^M!l(0))" GL-iCtlli^l 

(29) 



^^v ' / V ' / r/ 



where we have implicitly assumed that the boundary matrix Q_^ is non singular. 

(- s ■) 
By substituting the expression for the weights £ A' ^ (A) in the above equation and after 



several manipulations we find that Ti can be written as 

L-2 3 



« = ^EE 



7] 



i=l a=l 

2 3 r 



1 / 1 



1)2 



ri\S+\ 



2 3 






O 



L-1 ']Sl^Sl\ ^L+l 



V \s+l I ^ 



Q = l 



crUSl + S2^U + 2al_,{Slfal^, + ( _ _ 5(5 + 1) ) al_,al^. 



1 



(30) 



where the index a = 1,2,3 means the angular momenta components x,y,z, respectively. In 
the expression (jHUjl is also assumed the following toroidal boundary condition between the sites 
1 and L + 1 mi 



/ 



a 



L+l 



\ 



Oj 






/ 



1 

fi'llfi'22 — 9\2921 



all 



-all 



\ 



-9I2 



9I2 



/.+ \ 



a 



(J, 



(31) 



-9ii92i 

912922 

\ -'^911912 '2921922 911922 + 912921 I 

We see that the boundary matrix ^^ permits a more general type of interactions between 
the bulk spin-1/2 and the spin-S" impurity. The model 



\^l J 



turns out to be an interesting 
extension of the impurity system originally proposed in ref. jT2]. Finally, the eigenvalues E'^^'^ = 



f,\n{Alt-'^-'\x) 



|A=0 



of this Hamiltonian are given by 



E(±) 



1 
V 



L-1 



S + 



1 "^ 



(±)i2 _ af. 
4 



r[A, 



(32) 



where A^ satisfy the same Bethe ansatz equation 



2.2 The non-compact SU (1,1) realization 

Another possible realization of the Yang-Baxter algebra (Jl for the SU(2) i?-matrix is in 
terms of the SU{1, 1) Lie algebra. One motivation to study integrable models whose quantum 
space are SU{1, 1) invariant comes from the physics of atoms interacting with electromagnetic 
fields. In fact, different realizations of this symmetry in terms of bosonic operators lead us 
to a number of solvable atom-fields models ^3 1^ that are interesting generalizations of the 
famous Jaynes and Cummings paradigm ITT]. In this context, the study of these systems with 



general boundary representations appears to be important, since recently it has been argued 
[TH] that suitable combinations between boundary matrices ^^ and /^-operators can generated 
the so-called counter-rotating terms which are relevant in the case of high intensity fields. 
The /^-operator realization in terms of the SU{1, 1) generators K^ and K^ is given by J19j . 



where kj is the Bargmann index |20] which characterize the unitary representations of the 
SU{1, 1) algebra and c is any complex constant. The commuting transfer matrix of the corre- 
sponding mixed vertex model is therefore, 

T^%'-''^\X) = Tr^ [gAC±''^ . . . 4i'''^] • (34) 

The spectrum of this operator will depend much on the kind of representation we choose for 
the quantum space. One relevant representation for our purposes is the positive discrete series 
D~^{kj) where kj assume integer or half-integer values. More specifically, let the ket \kj,nj), 
with Uj = 0, 1, 2, . . . be the basis of D^{kj), then the action of the generators is given by 

K-\kj,nj) = {nj + kj)\kj,nj), (35a) 

K+\kj,n,) = ^{n, + l){n,+2k,)\kj,nj + l), (35b) 

K^ \kj, rij) = \/nj{nj + 2kj — 1) \kj, Uj — 1) . (35c) 



As before, we should now look for a local reference state \0)\ E SU{1, 1) such that the 

—1 (- k) 
ion of the transformed Mj'^-^^CJj ^ (A)Mj on it gives at 

This condition implies the following annihilation property 



—1 (- k) 
action of the transformed Mj'^-^^CJj ^ (A)Mj on it gives as a result an up triangular matrix. 



{[A + c] (p,+i - p,) - vK]ip,+i + pj) + vK^Pj+iP, + TiKj] |0)f ^ = 0. (36) 

In order to solve Eq.()36|) we shall repeat the procedure of previous section, i.e. we expand 

oo 

the pseudovacuum \Q)- = 2, C{kj,nj) \kj,nj) in the SU{1, 1) basis and use the properties 
()35all35c|) as well as that pj is the same for all j to generate relations for the coefficients 

10 



C{kj,nj). These relations are disentangled recursively and we find that the appropriate state 
|0>f'' is 



and that the action of the ^-operator on it is given by 

Before proceeding it is important to remark that the norm of \0)j is \C{kj, 0)^(1— p^)"^'^^ 
which means that this state can be normalized except when pj = ±1. This provides an extra 
requirement to the entries of the boundary matrix due to the constraint (j^n|) . Apart from 
this fact the next steps to determine the eigenvalues and eigenf unctions of T^i'"' ^ i^) by the 
quantum inverse scattering method are fairly parallel to those already described in section ITTl 
Omitting these details the final result for the eigenvalues is 

i=l i=l \ - A - c 



+^iji n(^ + - + ^^.o n '/\(.!'^^ (39) 



provided that the rapidities satisfy the Bethe ansatz equations 



L 



'>^'-vkA 9iJ^«A<«-Ar+„ 



where we have performed the convenient shift A^- — ^ A^^ — c. 

We would like to finished this section with the following comment. Having at hand two 

different families of /^-operators £)^j(A) G 1 I ®V^ k = 1,2 which satisfy the Yang-Baxter 

algebra with the same i?-matrix, the co-multiplication structure of this algebra allows us to 
construct even more generalized mixed vertex models by combining the tensor product of these 
two possible realizations. The corresponding monodromy matrix can be written as 

r^iX) = Ga^al,+l,{X)Cal,+l,--i{X) ■ . .Z^2(A)£^i(A), (41) 

11 



where CyXj{\) is defined by 



4'!(A), ifjG{7i,...,7Lj 



^aM) = { ^^l ' "^^ (42) 

£j4J(A), otherwise, 

and the partition {71, . . . ,7i} denotes a set of integer indices assuming values in the interval 

1 < 7i < ^1 + U- 

For example, one choose as the first family of /^-operators the higher spin operators Q and 
as the second one the 5'[/(l, 1) realization (jHH|) . Clearly, the eigenvalues problem associated 
with the general monodromy matrix (PT|) is solvable by trivial combination of the results of 
section 121 The global reference state |0) turns out to be^ 

io)(^)= n ®io)f^ n ®io)f^ (43) 

jG{/3i,...,/3i} i . 

otherwise 

while the eigenvalues are 

Ai - A + r/ 



A({5.},{..})(A) = ,w n (a + 1 + r^5.) n (^ + ^ - ^^^o n 



j=l j=l i=l '^i ^ 

i=i i=i i=i ^ '^^ 



+4^ n (^ + i - ^^.•) n (^ + - + ^^^) n \_\(.t' ' (^4) 

and Aj satisfies the Bethe ansatz equation 






3 SU{N) mixed vertex model 



A natural extension of previous section would be to consider mixed vertex models whose 
Boltzmann weights that combines isomorphic and non-isomorphic /^-operators based on the 



^We recall that here we are supposing that the boundary elements respect the condition pj 7^ ±1 to make 
the state |0)j' normalizable. 
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SU{N) symmetry at different orders of representations. Previous results for this mixed model 
[21], however, suggest us that the details entering the solution of such systems with general 
twists will follow closely that already described in section ITTl and in our previous work |llj . 

However, because SU{N) is not self-conjugate for A^ > 3 one expects that the Yang- 
Baxter algebra (^ admits further classes of realization on this group other than the higher 
spin representations. The purpose of this section is to explore such possibility by considering 
>C-operators whose quantum space is invariant by the conjugate representation of SU{N). 
Originally this /^-operator was discovered in the context of factorizable theories, representing 
the scattering matrix between particles and anti-particles ^2|- Recently ^^, this amplitude 
has been discussed more generally in terms of the braid-monoid algebra. We also recall that 
mixed vertex model combining the fundamental and conjugate representations are of direct 
interest in statistical mechanics since it appears to be related with the combinatorial problem 
of coloring the edges of the square lattice and fully packed loop models J2l]. The commuting 
transfer matrix T^^^'^^\\) of such mixed vertex model with a boundary Q_a is defined by 

Eqs. fj5l7p . namely 



Eqs. fj3l41|) . The first CA{\) operator is the fundamental SU{N) realization obtained from 



N N 



^%w = ^ E ^- ® ^- + ^ E ^S ® 41 (46) 

a=l a,/9=l 

while the second one CA(X) intertwines between the fundamental and conjugate representation 
of SU{N) and its expression is [22] 






N N 

^a,l3 ® ^N+l-a,N+l-l3- i'^V 

a,/3=l 

where p is an extra free parameter. 

Our first task is to find the reference states for each transformed operator £^(A) = 
M~^^Cj^-{X)Mj that annihilates all the ^ ^~ ' lower left elements for arbitrary A. The form 
of such pseudovacuums has a direct dependence on the elements of the gauge matrix which 

N 

here will be denoted by Mj = y, "^j(q^)/5)^o/3 • The first local reference state \0)- has been 

a,f3=l 



13 



already determined in ref. JT] and it is given by the first column of the gauge matrix Mj, 
namely 



( m,(l,l) ^ 
m,(2,l) 



10)^ 



(48) 



mj[N - 1,1) 
provided that the following ratios are satisfied 



Po^fi 



for a,/3 = 1,...,A^-1. 



(49) 



m,-(iV,/3) mj+i{N,f3y 

The second reference state |0) J. has a more involved representation in terms of the elements 
of My It turns out that this reference state can be put in the following form 



|0)f 



^N-1,N 



Si) 



(50) 



\ l,A^ / 

where c^l^ are the cofactors of the gauge matrix Mj. These are simply obtained from the (a, (3) 
minors of Mj with an appropriate sign, 






^a+|3 



mj(l,l) 



mj{a — 1, 1) 
mj{a + 1, 1) 

m,(iV,l) 



m,(l,/5-l) mj(l,/3 + l) 

mj(a — l,/9— 1) mj(a — 1,/? + 1) 

mj(a + l,/3- 1) mj(a + l,/? + l) 

mj{N,[3-l) mj{N,(3 + l) 



. m,(l,iV) 

. mj{a — 1, A^) 
. mj{a + l,N) 



mj{N,N) 



(51) 
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It follows that the action of the transformed /^-operator for fundamental and conjugate 
representation on their respective reference state is given in terms of an up triangular matrix, 



( a?(A)^ 



>cS(A) |o)5 



v 



where /^ are given by 








# 








# 
# 



# 



\ 



1 




aAr_l(A)7FfT- # 



|0)?, (52) 






V = < 



mj{N,a), a = l,...,N — l 



\y[rn,{N,z)^ 



det(Mj), a = N, 



(53) 



W, 



while the rational functions % (X) are 



4'Ha) 



A + ?7 for k = 1 
X for A; = 2, . . . , AT 



,(2) 



ar(A) 



^ - p for A; = 1, . . . , iV - 1 
^ - p - 1 for A; = A^. 



(54) 



The next step is to use the the remaining freedom of the elements of the gauge matrix to 
transform M^^Qj^M^^i into a diagonal matrix. This condition together with Eg . (1491) impose 
severe restrictions on the possible values for the ratios p^,/? which turns out to be the same 
satisfied by the ratio of the components of the eigenvectors of the boundary matrix Q^. This 
results in N possible choices for p^.^, / = 1, . . . , A^ and therefore from Eqs. ()48f5()j) there exists 
A^ type of appropriate local references states that will be denoted by \0y^' \ As a consequence 
of that we have A^ possible choices for the global pseudovacuum which can be written as 



|0) 



(0 



n ®io)r^ n ®io)f''^ for /=i 



AT. 



(55) 



i6{/3i,...,/3i,} 



otherwise 

Further progresses are made by choosing a suitable representation for the gauge transformed 
monodromy matrix that are able to distinguish the creation and annihilation fields over the 



15 



reference state fl55|) . Previous experience with vertex models having the triangular property 
(jHH) suggest us to take the structure used in the nested Bethe ansatz formulation of periodic 
SU{N) models 01111! 



/ 



X4(A) 



A{X) 5i(A) 

Ci{X) Du{X) 



Bn-i{X) 

DlN-l{\) 



\ 



\ C^-iiX) DN-niX) 



D 



N-lN-l 



(56) 



W J 



NxN 



By comparing Eq.()52|) with Eq. lj^UI) we observe that the operators Bk{X) produce new 
states when applied to the pseudovacuum state |0) while Ck{X) are clearly annihilation fields. 
Furthermore, we can also read the action of the diagonal fields, 



A(A)|0)« = [a?{X)]'^^[af\X)f^^\0f\ 

5..(A)|0)« = [aSi(A)]"M«Si(A)]^^^|0)(') 

/fc+i 



k = l,...,N-l. 



(57) 
(58) 



We now seek for further eigenstates |0) of T^^'^'^^\X) with the following structure 



ao 



,(1,0 



{l,/)x^'^„(0-'^l ,„,(/) 



5.,(Ar)...i?a (JAlTo)-^'"^ |0) 



m['^ ' m\ 



(59) 



a (i)...ai 

where the indices aj of the coefficients J-' ""i run over A^ — 1 possible values. The rapidities 
A,- will be determined by solving the eigenvalue equation 



T(^^'^^)(A)|0)(')=A(^^'^^)(A) 



dl) 



(60) 



At this point we have the basic ingredients to follow the general strategy of the nested Bethe 
ansatz method. We apply either A{X) or Dkk{X) on |</)) with the help of commutation rules 
that are the same known for the periodic case as explained in our previous work jTT]. From 
now on the computations are standard and since the few adaptations needed have already been 
described in ref.^J we shall present here only the final results. We find that the eigenvalues 
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^{L^M) 



(,) (A) are given by 



(LiM) 



,(i,0~ 



,(iv-i,0~ 



^L^o^" ™« (A; {A^ }, . . . , {Ar-^''^}) = 9^^^ (A + r^)- 

'"'1 ■•• '"']\j — i 



i.™S"A('.')_A + , 



" n^ 



i=i i 



A . -'^ - A 



'"^ A - Af '') + V ik >^f^ - A + r] 



■ ' ^ fc=i i=i ^ ~ ^j i=i ^'■- ~ ^ 



n ' 



'>n 



■"- A - Af -'•" + , 



r m 



(Af-1,0 ' 



(61) 



where the phase factors g^^'' are the eigenvalues of the boundary gauge matrix ^_4 ordered in 
such way that they satisfy the relation g^^' = g^^^^> for / = 1, . . . , A^. The corresponding nested 
Bethe ansatz equations can be written as 



,(0 



vC+l) 



A^'^+r^' 



A 



(1,0 



(0 



(0 



n 



■ Af'" - A"-" + „ '% Af -" - Aj-'" + „ 



n 



3 = 1 'j 



A^'^ - Af '') + r^ i\ \f'^ - Af •') 



(62) 



(0 



gii+k-i) "j-i yi^.i) _ xf-^.i) + ^ 



j(/+fe) 



i=i 






(0 
m) 
k 

n 



A<") - Af •'> + , 



™fc+l X (fc+1,/) _ X 



(fc,/) 



^=1 Aj- 



(fc,o 



A) ' + V 



n 



+ ^7 



A 



(k+i,i) .(fc,o 



A 



iV-2 



(63) 



9 



il+N-2) 



(l+N-l) 



9 



A 



K - vp 

(7V-1,0 



r]{p-l) 



(0 



-p-j- Aj - A,- + r] 



^N-1 ^iN-l,l) A(Af-l,/) 



i=i 



^(JV-1,/) _ ^^(N-2,l) 



n 

j=i 



A 



A 



J 



+ 77 



U.A^-1,0_^(A^-M)+^- 



(64) 



4 Quantum space transformations 

One striking result of previous sections is that the final forms of the transfer matrix eigenvalues 
and the Bethe ansatz equations are similar to that expected for vertex models with diagonal 
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boundaries. The corresponding diagonal twists are just the eigenvalues of the non-diagonal 
boundary matrix Qj^. This result motive us to look for possible correspondences between the 
eigenvectors of such vertex models with non-diagonal and diagonal boundaries as well. 

In order to establish the above mentioned relationship one has to explore the possibility 
of making quantum space transformations on the gauge transformed /^-operators. Rather 
remarkably, for all the vertex models discussed in this paper it is always possible to choose a 
invertible transformations Uj on the space Vj such that 

U-'CaAX)U, = CMX). (65) 

This means that quantum space transformations can be used to undo the modifications 
carried out on the /^-operators by the gauge transformations Mj. Let us now exemplify the 
importance of this property on the diagonalization of the standard transfer matrix defined by 
Eqs. (|3l6|) . Denoting by M^ the matrix that diagonalize the boundary matrix Q^ we can write 
T(A) as 

T(A) = Tr^[M^D^M_^iM^(M_^i/:^z.(A)M^)...(M//:^i(A)M^)M/], (66) 



Tr 



A 



D^Cj,l{\)Cal-i{\)...Cm{>^) , (67) 



where D^ is diagonal matrix whose entries are the eigenvalues of Qj\^ and the /^-operators are 
given by 

Ca,{X)=M^^CaAX)Mj,. (68) 

Now motivated by property ()f)5j) we can define a new operator T'(A) 

L L 

T'(A) = J] ®Ur^T{\) J] ®Uj = Tr^ [Dj,Cal{\) ■ ■ ■ C^W], (69) 

j=i j=i 

which is precisely the transfer matrix of the vertex model we have started with with diagonal 
twist Dj^. 

At this point is important to recall that the transfer matrix T'(A) can be diagonalized with 
very little difference from the periodic case because the diagonal boundary D^ does not change 
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drastically the properties of the monodromy matrix elements. This not only explain the reason 
why the Bethe ansatz equations and eigenvalues have the same shape for both diagonal and 
non-diagonal boundaries but also makes it possible to substantiate a clear relation between 
their eigenvectors. In fact, if we denote by \il)') one possible eigenstate of T'(A) then it follows 
directly from ()69|) that the corresponding eigenkets \'il)) of T(A) will be given by 

L 

and that the eigenvalues of T'(A) and T(A) are of course identical. 

Explicit expressions for the quantum matrices Uj are summarized in appendix B for all 
the vertex models discussed so far. It turns out that the above observations can be used to 
diagonalize the transfer matrix of even more complicated mixed vertex models. For instance, 
let us consider the following operator 



Tg{\) = Tr [^^(^)£^l(A)6?^(^-^)£^l-i(A) . . .Qa"^^^ Cai{\) 

that combines boundary matrices and /^-operators at any site of the lattice. 

Assuming that Qa are non-singular matrices one can insert Qa 
the trace ()71|). permitting us to rewrite the transfer matrix Tg{\) as 



Qa 



U) 



(71) 



terms all over 



T,(A) = Tr gA^'^^CAL{X)CAL-i{X) • • • i:^i(A) 



(72) 



where the new boundary matrix Qa is given in terms of the following ordered product 



Ga 



(ef) 



g/'^0/'-'^.--0A 



(1) 



(73) 



The transformed /^-operators in Eq. ()72|l are given by an extension of formula ()(J8|1 that 
accommodates site dependent transformations, namely 

-1 



CaAX) 



M 



A 



CaAX)m'X 



U) 



(i) 



in such way that for each site j the matrices M_^ are 



M 



(i) 



A 



.0--i)r7,(i-2) 



gA''~'>GA''''>---QA''>. 



(1) 



(74) 



(75) 
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As before, the next step is to find the appropriate quantum transformations Uj in order to 
reverse the gauge transformations M^ in Eq.()72|). Note that the matrix Uj are now different 
for each site j because its elements depend crucially on the entries of M^ , see appendix B for 
details. By performing the transformation (jU^ we are able to turn the transfer matrix Tg{X) 
into 

T;(A) = Tr^ [g/'^^CAL{\)CAL^i{\) ■ ..CAiiX)]. (76) 

Once again we end up with a transfer matrix problem with only one boundary matrix 
that has been discussed in details through this paper. Interesting enough, the diagonal twists 
entering into the Bethe ansatz equations and the transfer matrix eigenvalues are exactly the 
eigenvalues of the ordered product (ffH|) of all the boundary matrices. We remark that our 
findings when applicable to the two-site transfer matrix T^^' ^ i^) oi section 2.1 reproduce 
the eigenvalues and the Bethe ansatz equations proposed recently in ref . p[^ except for the fact 
that the number of roots vary up to 2(51 + 'S'2) instead of being fixed at this upper bound 
value. 

In what follows we will present two concrete examples of the utility of quantum space 
transformations in the problem of transfer matrices diagonalization. 



4.1 Spin-^ Heisenberg model 

The classical analogue of the integrable spin-S* Heisenberg model jHl E] is known to be a 25* + 1 
state vertex model whose Boltzmann weights are identified with the matrix elements of the 
following SU{2) invariant /^-operator [2113] 

25 

^S(^) = E/'(^)^'' (77) 

1=0 

and Pi is the projector onto SU(2)i in the Clebsch-Gordan 

decomposition SU{2)s ® SU{2)s- This operator is conveniently represented by the expression 

25 "* "* 
p^^^ S0S,-x, ^ (78) 

k=0 
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with Xk = llKl + 1) - 2S{S + 1)]. 

As usual, the transfer matrix associated to this SU{2) vertex model with an integrable 
boundary Qj( is 



(79) 



After some computation we find that the most general boundary condition compatible with 
integrability can be written as 



/ 



Ga 



9s,s 9s,s-i ■ ■ ■ 9s,-s 

gs-i,s gs-i,s-i ■ ■ ■ gs-i,-s 



\ 



^0) 



y 9-s,s 9-s,3-i ■ ■ ■ 9-s-s J 
where the coefficients gm,i satisfy the recurrence relation 

\/2MS^^m){S~+YTrri)gs-i,s9rn+i,i - \/2S{S - 1){S + I + l)9s,s-i9m,i+i 



9ni,l 



ioT I ^ m and l,m = —S, . . . , S and 



2S{1 - m)gs,, 



^1) 



9i,i 



for / 



9s~i,s^i9i+i,i+i \/2{2S -1){S -l-l){S + 1 + 2)gs-2,s9i+2,i + 2{S - I - l)gs-i,s9i+i,i 



9s,s 



-b, . . . , b . 



y/2b{b + l + l){b-l)gs, 



^2) 



As expected from the SU{2) symmetry this boundary matrix has four free parameters 
represented by the elements gs,s,9s,s-ij9s-i,s and gs-i,s-i- The standard procedure used 
to diagonalize transfer matrix based on higher spin representations consists in exploring its 
commutation with that of the mixed vertex model (fTUI) when 5*1 = 5*2 = ■ ■ ■ = 5^ = 5*. Unfor- 
tunately, this property is no longer valid for arbitrary non-diagonal boundaries and therefore 
another route has to be taken in these cases. One possible way is to proceed exactly as ex- 
plained in the beginning of section |3] We first find the matrix M^ that diagonalize the 
boundary matrix Q^ (jHOl) and afterwards we perform the quantum space transformation (jU3j) 
by using the matrices t/- collected in the appendix B. The transformed transfer matrix T^(A) 
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is then given by 



s ris)i 



T's{\) = Tt^ D'^C]^>{\)...C'^I{X) 



<S), 



13) 



where -D^ is the diagonal matrix having the following two possible forms 



n^ 

^^(H 



diagUW^r^:^ 



(9^ 



(-) 



S-m 



i9^^f-' 



{9. 



{-) 



2S 



(+) ' ■ ■ ■ ' 

9s 



(9. 



(+),5'-m-l ' • ■ ■ 



(9^ 



(+) 25-2- 



(9. 



(+),25-i ; ' 



M) 



and 



D 



■A{-) 



diag 



(9 



i-h'' 



(9 



(_) 2S-1 



(9 



fyS+m 



(9^ 



(-)^ 



A+h^^-^' 



(^) 2S-2- 



/_!_■) S+m—1 ' 



9?' 



.s'f\s'f 



») 



XgD (sD (si. 

The parameters (yf^. are two particular eigenvalues that are able to parameterize the total 
eigenspectrum of the boundary matrix ()80|1 . We recall that these two possibilities we have 
at our disposal is related to a remaining Z2 symmetry of the Hilbert space. Now Tg{X) can 
be diagonalized by the quantum inverse scattering method adapting the results of ref. j3] to 
include diagonal twists. It turns out that the final result for the eigenvalues of Ts{X) is 
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aS(a) = E 



^g(-))STm 



n± 



m=-S \9s ) 

where the functions tm(A) and g.m(A) are given by 

s 






(±) , V^ 



(86) 



tm{X) = {X + 2r]S) H 



fc=m+l 



\ + Tjk — TjS 

X + rjk + rjS' 



(87) 



and 



qm{X) 



{X + r]S + ^)iX-riS-^) 



(A + r]m + ^)(A + rjm — 2) ' 
while the corresponding Bethe ansatz equation for the rapidities A^ are 



,(±) 



XT' +js^ 

TjSj 



A 



(±) 






Ar' + „ 



a!*'-,/ 



») 
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4.2 Singular boundary matrix 

In the most part of this paper we have imphcitly supposed that the boundary matrices Qji 
are invertible. Here we would hke to explore the interesting situation in which these matrices 

become singular. For sake of simplicity we shall discuss this problem for the transfer matrix 

(i i) ___ (i i) 

^1/2 '^ (A) (fTUj) that are built up in terms of the simplest ^-operators jCJj^ (A). A possible 

representation of the corresponding singular boundary matrix is 

(„ 311922 

^" ^- I, (90) 

921 922 

where we are assuming that (721 7^ 0. 

For general values of Qij this matrix is still diagonalizable and it is similar to a diagonal 
matrix T>'^^^\ namely 

V^'3^ =1 I , (91) 

\ 5-11 + 5-22 / 

while the explicit expression for the matrix iVf4 that diagonalize Qj^^ is 

(322 511 
321 321 I _ (92) 

We now proceed by performing the gauge transformations ()67|) with the help of the matrix 
(|92j) as well as the quantum space transformation [\<6b\ where Uj = M^^ . As a result, the 
transformed transfer matrix we need to diagonalize becomes 



/(ov.o)/ 



)(S3) /'(2'2)/'\^ /'(2'2) 



T%-'--\\) = Tr^ [D'^'>C'J;^^'{X) . ..C]j;-^'{X)\. (93) 

Clearly, this transfer matrix is defective since T'^^s') contains one diagonal vanishing element 
and therefore it is proportional to a single diagonal monodromy matrix element. Thanks to 
this property, the non-null eigenstates \ip') of this operator can be found by direct inspection 
and they are given by 
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for n = 0, 1, . . . , L. The corresponding eigenvalues A^"-'(A) have the following form 

A(")(A) = (^11 + ^22)[A]"[A + r/](^-"). (95) 



eigenstates of T-^?2 "'^ (A) are 



Considering the result (jMj) and Eq. (f7n|) with Uj = M_^^ we conclude that the original 

2'' 
/2 

n I _922\ ^ I 911 \ 

(96) 




Finally, we remark that the results ()95l96p have been previously conjectured in ref. JT] on 
basis of exact diagonalization of T^J^''^ (A) for some values of L. This discussion not only 
clarifies the origin of such results but also emphasizes the importance of the quantum space 



transformations in the general framework of transfer matrices diagonalization. 

5 Conclusions 

In this paper we have investigated integrable non-diagonal toroidal boundary conditions for 
mixed vertex models and related spin chains based on the isotropic SU{N) i?-matrix. In 
particular, we have applied the quantum inverse scattering method to obtain the transfer 
matrix eigenvalues and the corresponding Bethe ansatz equations for systems with quantum 
spaces invariant by the high spin representation of SU{2), the discrete D~^{k) representation 
of SU{1, 1) and the conjugate representation of SU{N). 

We have introduced the notion of quantum space transformations that together with the 
standard Baxter's gauge transformations ^3] allowed us to diagonalize the transfer matrix of 
general families of vertex models that mix both boundary and /^-operator representations. We 
have used this approach to solve the Heisenberg model with arbitrary spin and non-diagonal 
twists as well as the case of singular boundary matrices. We strong believe that this method is 
easily applied to other isotropic vertex models such as those invariant by the 0{N) and Sp{2N) 
Lie algebras |2S]. We also hope that this property could be generalized to accommodate the 
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solution of trigonometric vertex models and therefore improving our understanding of the 
integrability. 
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Appendix A: Completeness for SU{2) mixed vertex mod- 
els 



In this appendix we consider the completeness of the Bethe ansatz solution of the SU{2) 
invariant vertex models discussed in section ITTl for L = 2. More precisely, we are going to 
verify that all the (25*1 + 1)(2S'2 + 1) eigenvalues of the transfer matrix T^i2 W P^ can be 
obtained from the Eqs.(I^ 



for some values of 5*1 and 5*2: 



.(ii), 



1/2 (A) 

We can search for solutions of the Bethe ansatz equations ()28|) beginning either with |0) 
or |0) eigenvectors. By choosing the state |0) the zero-particle eigenvalue can directly be 
read from (j27j), 

A?''\a) = ^;;)(A + r^)(A + ^) + ^;7)(A)(A - |). (A.l) 

For the one-particle state B{\[ ) |0) we find two possible Bethe ansatz roots given by 






and the respective eigenvalues are 

-.2' 



aSI'^V) 



yi/2 



A^ + ?A + 



4 



^yi/2 



4 U]i 



y + -ix 



yi/2 



(A.2) 






± 



(^i72+^i;2^)^+ 



32g[-}lg[-l (A.3) 
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The next step is to consider the two-particle state B{\\ )B{X2 ) |0) and again two 
possibihties are obtained and they are 



A 



(+) 

i± 



5^(^1/2 



J+) 

yi/2 



yi/2 yi/2 



T 



8(.S-.!;20 



r/v^^5((?;/ - (?W') - 38^;-)^;;) ± 3 (^g[-l + g[y,) ^{g^^l + (?J+V + 32^^^ 



1/2 



yi/2 



yi/2 



(A.4) 



and 



A 



{+) 

2± 



5^ (^V2 + 9yl 

8 Ul/2 - 9i/2 



T 



iaijl 



& 



+ ^^&yl 



sig 



'1/2 



J+) 

yi/2 



rj^/2^5{g[jf - g[f) - 38g^^lg[-}l ± 3 (^g[-^l + g[}l) ^{g[-^l + g\)lf + 32^i7)^i; 
The corresponding eigenvalues are then given by 



(A.5) 



A, 



(ii)^^- J-) 



2± 



\>^)=9V/2 



A^ + ?A + ^ 



yi/2 



A^ + ?A-^ 



±\/ (^1/2 +^1/2 



(+)a2 



32(?;+^(?jT,\ (A.6) 
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The last sector consists of three Bethe ansatz roots and an analytical expression for them 
has been prevented within our numerical resources. We have however performed an extensive 
numerical work confirming that the associated eigenvalue is 



T, 



(1,1), 



A?''\a) = g^-/l{\ + r])(A + 1^) + ^W(A)(A - |). 

1/2 '(A) 
As before, starting with the state |0) it is not difficult to see that 

The next state is the one-particle sector giving us two possibilities 



(A.7) 



(A.8) 







'AT. ± 



'U 



9l/2 =F A/fi'i/2 



(A.9) 



26 



and the corresponding eigenvalues are 



Ai^i^^A) = ,W(A + f )(A + |) + 4^(A - |)(A + f ) ± 4^)1^,. 



(A.IO) 



On the other hand the two-particle state provides us three types of Bethe ansatz roots. 

The first one is 

S-) ^ J+h 



\(+) (5'i/2 +fi'i/2) 1 x(+) n 

Ai = V-jT) ^ and A^ ' = 0, 



(^172-^172) 



and the respective eigenvalue is 



Ar'(A) = (9:>9!7.')(^ + f)(^-i)- 



The other remaining solutions are given by 









TV 






+r] 



1^1/2 



yi/2) 






^1/2 



1^1/2 



^«)2 + 32^W^i-) 



2v^(^l72 



yi/2^ 



and 



^2± 



%i72-^i7J) ^"^ %!72-^i;j) 

^(^172 + ^172)' - 1<2^172 ± (^172' - 9[f)^{9[-}l + 9[/ly + ^2g[yA/l 



2j 



2V2{g[/i - (?;; 



while the associated eigenvalues are 



1/2 



(A + ^)(A-^) + ^ 



± 



(^2 +^172)^ + 32^^72 • 



(A.ll) 



(A.12) 



(A.f3) 



(A.14) 



(A.f5) 



The remaining sectors have been investigated numerically and we verified that the three 
and the four particle states produce the following eigenvalues 



M^)|^^-J-)|^ , 3^^M , ^a , J+)m ^^M , ^^ ^ . /.(+).(-) 



A;,r (A) = ^172^(A + f )(A + ^) + ^i72^(A - ^)(A + ^) ± y'^i;^^!/ 



/2' 



(A.16) 
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and 



^12 



'(A)=^K[A+^r+^i;iA-^r- (A.17) 

Finally, we remark that as byproduct of this analysis we have also checked the completeness 
of the eigenspectrum of the transfer matrix T| ' (A) of section 14.11 since its Bethe ansatz 
equation are the same as that discussed above. 



Appendix B: Quantum space transformation 

In this appendix we will present the explicit expressions for the quantum space transformations 
Uj that are able to undo the gauge transformations on the auxiliary space. The matrix elements 
of Uj are directly related to entries of the gauge matrix M^ since they have been required to 
satisfy the following identity 



CA,{\) = u-'[M^iY'^A,WM^i^a 



V 



(B.i: 



for each /^-operator mentioned in the main text. In order to present this relationship it is 
convenient to write the gauge matrix as 



N 



M 



U) 



A 



a,/3=l 



- (i) "{A) 



(B.2) 



where we recall that A^ is the dimension of the auxiliary space A. We now start to list the 
quantum matrices Uj for the /^-operators used in this paper 



YotC%'^\x): 



( 



U, 



(i-s,) 



^s„s, 



Us,,S,-l 



Sj l,Sj Sj l,Sj 1 



where some of the matrix elements Uij are 



s,-i 






(B.3) 



Us,,S,-k = Ri.i 



;^{j)]2S,~k\-{j)-{k 



m 



(25, 



1,2J 



A;!(2S',-A;)!^ 



(B.4) 
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while the remaining ones satisfy the following recurrence relation 



.(i)™(i) ™(i)™(i) 



Us,-k-l,S, 



{mY{m^2^^ - mi^^m^^i) / n{2Sj - n + 1) 



^Uh 



\m 



i,iJ 



{k + l){2Sj-k) 



US,-k,S.i~n+l 



^(i)' 



+ 



m 



2,1 



k{2Sj -k + 1) m^il 2{Sj - k) 

Us,~k+l,Sj~n + 



^(ij; ]J (k + l){2S,-k) 
and k, n are integers satisfying k,n = 0, . . . , 2Sj. 



rr^W ik + mS.-k) '''^-'^'^- ^^-'^ 



For Cp\\y. 



u. 



where the first row are given by 



Uo,n 



(l.fc.) 



' Uofi Uo,l Mo,2 

Ul,0 Ul,l Ul,2 

U2,0 U2,l U2,2 

y : ; ; 



\ 



^or 



-D" 



m 



2,1 



(2% + n-l)! 



^l V n!(2A;,-l)! 



-^0,0 



m 



2,2 



and the other elements are obtained recursively by the following relation 



^(i)^(i) ^(i)^(i) 



Wn+l,; 



{mYim^^!^ - m^l^m)^}) j l{2kj + 1-1) 

{n+ l){2kj +n) 



(i)i2 

2,2J 



Un,l~l 



m^\ / n(2A;j- + n - 1) m^ 2(A;j + n) 



UnA, 



and n, / are integers n, I = 0,1,2, ... . 



(B.6) 



(B.7) 



fB.8l 



(s). 



For the higher spin operator Cj{X): 



Uf = Mf 



(B.9) 



(1)/ 



For the fundamental SU{N) operator £^.(A): 



t/j^) = M^i^ 
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For the conjugated SU{N) operator C%{X): 



( T^U) 



u. 



(2) 



.i3) 



^N,N ^N,N-l 



M) 



-Sj) 



'^N-1,N '^N-1,N~1 



\ H,N '^1,N-1 



where c^^L are the following cofactors 



^(j) 






\Q+/3 



m 



1,1 



m 



m 



U) 

(J) 
a+1,1 



^(j) 



m 



NA 



^U) 



m 



1,(3-1 



M) 



/3) \ 

^N,l 

.i3) 



'-1,1 / 



^ii) 



m 



l,/3+l 



M 



"^a-1,/3-1 "^a-1,/3+1 

- (i) - (i) 

''Wl,/3-l ''Wl,/3+l 



^(i) 



^(i) 



"^Af,/3-l "^Af,/3+l 



;^0-) 



m 



m 



1,N 



a-l,N 



m 



a+l,N 



^0) 



m 



N,N 



(B.ll) 
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